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Abstract. A unified first law of black-hole dynamics and relativistic thermodynamics 
is derived in spherically symmetric general relativity. This equation expresses the gradi- 
ent of the active gravitational energy E according to the Einstein equation, divided into 
energy-supply and work terms. Projecting the equation along the flow of thermodynamic 
matter and along the trapping horizon of a black hole yield, respectively, first laws of rel- 
ativistic thermodynamics and black-hole dynamics. In the black-hole case, this first law 
has the same form as the first law of black-hole statics, with static perturbations replaced 
by the derivative along the horizon. In particular, there is the expected term involving 
the area and surface gravity, where the dynamic surface gravity is defined by substitut- 
ing the Kodama vector and trapping horizon for the Killing vector and Killing horizon in 
the standard definition of static surface gravity. The remaining work term is consistent 
with, for instance, electromagnetic work in special relativity. The dynamic surface gravity 
vanishes for degenerate trapping horizons and satisfies certain inequalities involving the 
area and energy which have the same form as for stationary black holes. Turning to the 
thermodynamic case, the quasi-local first law has the same form, apart from a relativistic 
factor, as the classical first law of thermodynamics, involving heat supply and hydrody- 
namic work, but with E replacing the internal energy. Expanding E in the Newtonian 
limit shows that it incorporates the Newtonian mass, kinetic energy, gravitational poten- 
tial energy and thermal energy (internal energy with fixed zero). There is also a weak type 
of unified zeroth law: a Gibbs-like definition of thermal equilibrium requires constancy 
of an effective temperature, generalising the Tolman condition and the particular case of 
Hawking radiation, while gravithermal equilibrium further requires constancy of surface 
gravity. Finally, it is suggested that the energy operator of spherically symmetric quantum 
gravity is determined by the Kodama vector, which encodes a dynamic time related to E. 
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I. Introduction 

Analogies between the laws of thermodynamics and certain properties of black holes have 
aroused great interest, particularly since Hawking showed that the analogy becomes actual 
in at least one respect, namely that the surface gravity determines the effective temperature 
of a quantum field on a stationary black-hole background. This unexpected link between 
two hitherto unrelated branches of physics has inspired the hope that they may be unified. 
It is generally thought that there should be some deeper explanation involving a quantum 
theory of gravity. However, this work reveals a purely classical connection. 

Understanding this requires a paradigm shift from black-hole statics to black-hole dy- 
namics, meaning a theory of dynamically evolving black holes rather than just stationary 
black holes (and perturbations thereof). This is analogous to the shift from thermostatics 
to thermodynamics, meaning the theory of dynamically evolving thermal systems rather 
than just those in thermal equilibrium (and perturbations thereof). Many researchers ap- 
parently still believe that the relevant generalisation concerns event horizons, following 
the textbook definition of black holes by event horizons. However, this definition makes 
essential use of the global assumption of asymptotic flatness. The real universe is thought 
not to be asymptotically flat, so that event horizons do not actually exist. Moreover, even 
theoretically, the global nature of event horizons means that they cannot be physically lo- 
cated by observers. Their location is unknown until the universe has ended. For the same 
reason, event horizons have no local dynamical significance. An event horizon might be 
passing through the reader right now, without any physical effect. In short, it seems that 
event horizons are not the relevant generalisation. According to this view, the relevant 
second law of black-hole dynamics cannot be what the textbooks state, namely Hawking's 
result that the area of an event horizon is non-decreasing. This is reflected by the ab- 
sence of a corresponding first law for event horizons. The textbook first and zeroth laws 
concern stationary black holes, so are more accurately described as laws of black-hole stat- 
ics. The textbook first and second laws do not even involve the same derivative, thereby 
undermining the analogy with thermodynamics. 

A new framework for black holes was introduced in a previous paper [1], referred to 
as Paper I. This dispensed with global assumptions and proposed a general, dynamical, 
quasi-local definition of black holes in terms of trapping horizons, hypersurfaces foliated 
by marginal surfaces. Paper I also derived a second law of black-hole dynamics, expressing 
the increase of area element along the trapping horizon. This should not be confused with 
Hawking's second law for event horizons. In the new framework, black-hole dynamics is 
formulated as the dynamics of trapping horizons. 

A corresponding first law of black- hole dynamics requires a definition of surface gravity 
for non-stationary black holes. This article gives a new definition of surface gravity for 
spherically symmetric black holes, with the same form as the usual stationary definition 
involving the Killing vector on the Killing horizon, but instead using the Kodama vector 
on the trapping horizon. The Kodama vector encodes a preferred time in spherically 
symmetric space-times and has physical properties analogous to those of Killing vectors. 

Another key concept is energy. In general, there is no agreed definition of gravitational 
energy in the theory of general relativity. However, in spherical symmetry there is the 
Misner-Sharp energy E, which does have all the physical properties one might expect of 
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active gravitational energy, as explained in another recent paper [2] , referred to as Paper II. 
In particular, E has the expected behaviour in various limits: vacuum, small-sphere, large- 
sphere, Newtonian, test-particle and special-relativistic limits. It is also intimately related 
to the characteristic strong-gravity phenomena, namely black holes and singularities. 

The first law of black-hole dynamics given here just expresses the derivative of E along 
the trapping horizon, according to the Einstein equation. This turns out to involve the area 
and surface gravity as expected, with additional terms depending on the matter model, 
such as an electromagnetic term which has the form of electromagnetic work. This first 
law has the same form as the first law of black-hole statics, with the perturbation being 
replaced by the derivative along the horizon. In comparison to most references, it should 
be noted that this involves the energy on the horizon, rather than the energy at infinity. 
The new first and second laws involve the same derivative, the dynamical derivative along 
the trapping horizon. 

The next observation is that, for a perfect fluid, the derivative of E along the fluid 
flow is given by the hydrodynamic work, as shown in Paper II. This therefore reads as a 
relativistic version of the classical first law of thermodynamics, albeit for the thermody- 
namically trivial case with no thermal flux. To include such thermodynamic effects, it is 
necessary to generalise the model of a perfect fluid. A generally relativistic theory of ther- 
modynamics has been developed in an accompanying paper [3], referred to as Paper III. 
This theory models heat relativistically by an energy tensor, called the thermal energy 
tensor, whose various components are thermal energy density, flux and stress. (Thermal 
energy is internal energy with fixed zero) . This turns out to add a heat-supply term to the 
first law, consistent with the classical first law of thermodynamics. Actually, this requires 
no thermodynamics other than the form of the energy tensor given by Eckart [4] and the 
definition of heat supply as an integral of thermal flux. Thus we have a quasi-local first law 
of relativistic thermodynamics. This shows that the gravitational energy E incorporates 
the thermal energy, as may be verified explicitly in the Newtonian limit. Then the unified 
first law is just the equation expressing the gradient of E according to the Einstein equa- 
tion. Projecting the unified first law along the material flow or along a trapping horizon 
yields the first law of thermodynamics or black-hole dynamics, respectively. 

One might also expect there to be a unified zeroth law, concerning gravithermal equi- 
librium. The Gibbs definition of thermal equilibrium is constancy of the partial derivative 
of entropy with respect to thermal energy. Replacing thermal energy with gravitational 
energy E in this definition yields constancy of an effective temperature, differing from the 
actual temperature by a red-shift factor. This generalises the Tolman condition [6] for 
gravithermal equilibrium and agrees with the case of Hawking radiation [5] on a stationary 
background. Surface gravity is also constant for stationary black holes, thereby linking the 
two zeroth laws. 

The article is organised as follows. Section II reviews the relevant properties of spher- 
ically symmetric space-times, defining E and deriving the unified first law. Section III 
considers work in the case of electromagnetism and Section IV reviews black-hole stat- 
ics. Section V describes the Kodama vector or dynamic time, Section VI introduces the 
dynamic surface gravity and Section VII derives the first law of black-hole dynamics. Sec- 
tion VIII derives the corresponding first law of relativistic thermodynamics and shows 
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how E incorporates thermal energy. Section IX considers thermal equilibrium, gravither- 
mal equilibrium and the zeroth laws. Section X concludes. 

II. Unified first law 

In spherical symmetry, the area A of the spheres of symmetry is a geometrical invariant. 
It is convenient to use the areal radius r = a/ ' A/At:, giving 

A = 4nr 2 . (2.1) 

Following Papers I and II, a sphere is said to be untrapped, marginal or trapped as V"r 
is spatial, null or temporal respectively, V"r = being degenerate. Here V denotes the 
covariant derivative operator and denotes the contravariant dual (index raising) with 
respect to the space-time metric g. Similarly, b will denote the covariant dual (index 
lowering). If the space-time is time-orientable and V*r is future (respectively past) causal, 
then the sphere is said to be future (respectively past) trapped or marginal. Untrapped 
or marginal spheres have a spatial orientation given by the direction of V"r: a spatial or 
null direction is said to be outward (respectively inward) if r is increasing (respectively 
decreasing) in that direction. A hypersurface foliated by marginal spheres is called a 
trapping horizon. A trapping horizon is said to be outer, degenerate or inner as V 2 r > 0, 
V 2 r = or V 2 r < respectively. 

Paper II proposed that future (respectively past) outer trapping horizons be taken as 
the dynamical definition of the outer boundaries of black (respectively white) holes. Black 
holes are thereby characterised quasi- locally, actually locally in spherical symmetry. This 
definition should not be confused with those of event horizons or apparent horizons [7], 
both of which require the global assumption of asymptotic flatness. It is also unnecessary 
to make assumptions about the global nature of the trapping horizon, for instance that 
it remains of the future outer type. Such questions are still of interest, but depend on 
the dynamics of the chosen matter model. The results of this paper are essentially local, 
independent of the matter model except for energy conditions and, with minor exceptions 
for degenerate cases, apply to any trapping horizon. 

The Misner-Sharp energy [8] may be defined by 

E = \r (1 - Vr • V H r) (2.2) 

where the dot denotes contraction, the sign convention is that spatial metrics are positive 
definite and the Newtonian gravitational constant is unity. Paper II investigated the 
physical properties of E in detail, establishing that it represents active gravitational energy. 
This point deserves emphasis, since the literature contains many definitions of so-called 
energy which do not have the relevant properties in physically understood limits; a list of 
references may be found in [9] . In comparison, E does behave as active gravitational energy 
in the vacuum, small-sphere, large-sphere, Newtonian, test-particle and special-relativistic 
limits. 

Two invariants of the (contravariant) energy tensor T will be useful: the function 

w = -±traceT (2.3) 
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and the vector 



ij) = T ■ Vr + wVV 



(2.4) 



where the trace refers to the two-dimensional space normal to the spheres of symmetry. It 
is also convenient to use the areal volume 

V = fyrr 3 (2.5) 

due to the relation 

W = AVr. (2.6) 
Then the gradient of E is determined by the Einstein equations as 

VE = A^ + wVV. (2.7) 

This equation turns out to be the unified first law, as explained subsequently. This may 
be generalised beyond spherical symmetry in terms of the Hawking energy [10]. 

One may physically interpret w as an energy density and ip as an energy flux or 
momentum density. Actually, assuming the null energy condition, w is a lower bound for 
the energy density measured by an observer. Assuming the dominant energy condition, 
w > 0. Assuming the null energy condition, if) is past (respectively future) causal in future 
(respectively past) trapped regions, and outward spatial or null in untrapped regions. If 
the space-time is asymptotically flat, one finds that ip becomes future (respectively past) 
null at future (respectively past) null infinity, reducing to the Bondi flux. This can be seen 
from the first law (2.7), since wW tends to zero at infinity. Thus ip is a local version 
of the Bondi flux, actually the outward flux minus the inward flux. Then the first term 
Aip in the unified first law may be interpreted as an energy-supply term, analogous to the 
heat-supply term in the classical first law of thermodynamics, while the second term wW 
may be interpreted as a work term. 

To check the above properties, one may use double- null coordinates (£ + , £ _ ), in terms 
of which the line-element is 



ds 2 = r 2 d0 2 + 2# + _d£+d£" (2.8) 

where dfi 2 refers to the unit sphere and r > and g_| < are functions of Then the 

coordinate forms of (2.2-4) and (2.7) are 

E = r{\- g + ~d + rd-r) (2.9a) 

w = -g+-T+- (2.96) 

^ = T++d+rd + + T—d-rd- (2.9c) 

d±E = A^± + wd±V = -Ag+~ (T + _d±r - T ±± d T r) . (2.9d) 



The last expression was derived from the Einstein equations in Paper II. In terms of 
these coordinates, a sphere is trapped if d + rd-r > 0, being future trapped if d±r < 
and past trapped if d±r > 0, taking d± to be future pointing. In an untrapped region, 
d + rd-r < and one may locally choose the orientation d + r > 0, cLr < 0, meaning that 
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<9_i_ is outward and <9_ is inward. As in Papers I and II, note that the null energy condition 

requires T±± > and that the dominant energy condition requires > 0. The Einstein 

equations translated from Paper II are 

d±d±r - d± log(-g + -)d±r = -4nrT ± ± (2.10a) 
rd+d-r + d+rd-r - ±a+_ = 4yrr 2 T + _ (2.106) 
r 2 d+cL log(-o+_) - 2d+rd.r + g+_ = 8nr 2 (g+_T^ - T + _) (2.10c) 

where 9 is a standard latitude. 

III. Electromagnetic work 

Since the desired first law of black-hole dynamics will involve work terms, it is necessary 
to identify the correct definition of work for the relevant matter fields. This is done as 
follows for the Maxwell electromagnetic field (excluding magnetic monopoles). The only 
non-zero components of a spherically symmetric electromagnetic field tensor F are given 
by 

£ = -g + ~F + . = g+-F_ + (3.1) 

which may be intepreted as the electric field strength, the magnetic field vanishing. A 
standard definition of charge [11] is 



An 

where the integral is over a sphere of symmetry, the area form being implicit. Thus the 
field may be characterised by either S or e. The electromagnetic energy tensor T e has one 
independent component, yielding 

w e = S 2 /8n (3.3a) 
^ e = (3.36) 

where the notation is that of (2.3-4) with the subscript e referring to the electromagnetic 
field. Therefore the unified first law (2.7) reads 

VE = — + w VV + Aip b Q (3.4) 

where the subscript o refers to other (non-electromagnetic) fields, i.e. T = T e + T a and 
so on. In the absence of other fields, the Reissner-Nordstrom case is recovered. This 
identifies the electromagnetic work as S 2 W/8n. This agrees with the standard expression 
for electric work [12] in special relativity, which is 

W = — [ £ 2 dV (3.5) 
8?r Jy, 
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where dV is the volume form of a flat spatial hypersurface containing the region E. The 
integrand is just the local energy density w e of the electromagnetic field. Some authors 
accept this expression for work only for infinite E, allowing it to be rewritten in various 
ways involving boundary terms which vanish at infinity. However, such terms are generally 
non-zero at a regular centre. Indeed, typical suggestions become infinite. 

It is also possible to write the electromagnetic work in terms of the electric field 
covector 

£ = F-k (3.6) 

where k is the Kodama vector defined subsequently by (5.1). Then £ = £Vr, so that the 
electromagnetic work is 

— = - 2 e£. (3.7) 

The second expression has the same form as the work e$£ done on a test charge eo by 
the electric field £ in special relativity, with the half arising to avoid counting the self- 
interaction twice. In other words, the work is that done on the charge distribution by its 
own electric field. 



IV. First law of black-hole statics 

This section reviews the theory of black-hole statics, or equilibrium mechanics, which 
concerns stationary black holes; see e.g. Carter [13] or Wald [14]. Stationary space-times 
are defined by a Killing vector £ generating asymptotic time translations. That is, £ 
satisfies the Killing equation 

V®£ b = (4.1) 

where <g> denotes the symmetric tensor product. A Killing horizon is a null hypersurface 
such that £ is tangent to the null generators. Killing horizons define the boundaries of 
stationary black holes, white holes and cosmological regions. It follows from the definition 
that £ • (V A £ b ) is parallel to £ b on a Killing horizon, where A denotes the antisymmetric 
tensor product. The proportionality constant defines the surface gravity k: 

£ . (V A £ b ) = /cf b on a Killing horizon. (4.2) 



Taking the standard example of the source-free electromagnetic field, the only spherically 
symmetric solution to the Einstein-Maxwell equations is the Reissner-Nordstrom solution, 
parametrised by two constants, the charge e and asymptotic energy m; these provide an 
analogue of thermostatic state space for stationary black holes. For m > \e\ the solution 
describes a black hole with surface gravity k = r~ 2 \/m 2 — e 2 on the Killing horizon at 
radius r = m + \Jm 2 — e 2 . Evaluating these quantities on the horizon yields 

KdA ede , A „. 

dm = H 4.3) 

57T r 

where d denotes the differential in state space (m, e). This is the form usually given for the 
first law in this case, though it may be generalised to non-stationary perturbations [15] (of 
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stationary solutions). However, this form is not suitable for generalisation to space-times 
which are not asymptotically flat, due to the appearance of the asymptotic energy to. 
Instead, one may use the quasi-local energy E, which for the Reissner-Nordstrom solution 
is 

e 2 

E = m . (4.4) 

2r v ' 

Rewriting (4.3) in terms of E yields 

KdA + S 2 dV .... 
dE = (4.5) 

where 8 = e/r 2 as in (3.2). This is a physically acceptable form for the first law of 
black-hole statics in the electromagnetic case, since the last term may be interpreted as 
electromagnetic work, in agreement with special relativity (3.5) and the general expression 
(3.7) for electromagnetic work. Thus the analogue of internal energy for a black hole is the 
energy E on the horizon, rather than the energy at infinity; put another way, the energy of 
the black hole rather than the energy of the space-time. This is exactly as one would have 
expected by physical intuition and provides further evidence for the physical correctness 
of E as a definition of active gravitational energy. 

V. Dynamic time 

In black-hole dynamics, unlike black-hole statics, there is no Killing vector to define such 
quantities as surface gravity. However, in spherical symmetry there is a natural analogue, 
the vector 

k = curl r (5-1) 

where the curl refers to the two-dimensional space normal to the spheres of symmetry. In 
double-null coordinates, 

k = — g -1 (<9 + r<9_ — <9_r<9 + ) (5.2) 

up to orientation, which in an untrapped region may be locally fixed such that k is future- 
pointing. Kodama [16] introduced k and Paper II discussed some properties. It follows 
that 

k ■ Vr = (5.3a) 
l 2E 

k-k b = 1 (5.36) 

r 

which equivalently defines k. So k is spatial, null or temporal for trapped, marginal or 
untrapped spheres respectively. In particular, the definition of the boundary of stationary 
black holes, a hypersurface where the Killing vector is null, may be generalised to a hyper- 
surface where the Kodama vector is null, which is equivalent to the definition of trapping 
horizon. Unlike a Killing horizon, a trapping horizon is generally not null. 

The Kodama and Killing vectors agree for the Reissner-Nordstrom black hole. More 
generally, they agree in a stationary, spherically symmetric space-time if k commutes with 
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V*r. This condition implies that one may introduce coordinates (r, t) where k = d/dt, for 
which the line-element becomes 



ds 2 = r 2 d0 2 + ( 1 - —) \r 2 -(l-—\ dt 2 . (5.4) 



r J \ r J 

This is a Schwarzschild-like form, with E generally a function of (r, t). Then examine the 
various components of the Killing equation (4.1): 

i c d c g ab + 2g c(a d b) e = (5.5) 

for a Killing vector £ normal to the spheres of symmetry. The components tangential to 
the spheres of symmetry yield £ • Vr = 0, so that £ is parallel to k: £ = £°k. The rr 
component yields £ • VE = 0, the tt component yields k ■ V£° = 0, and the tr component 
yields V*r • V£° = 0, so that £° is constant. If the space-time is asymptotically flat and 
the Killing vector is normalised as usual by £ • £ b — > — 1 as r — > oo, then £° = 1 since 
• k 1 " — > — 1 as r — > oo. Thus the Kodama vector reduces to the Killing vector, £ = k. 
Another physically relevant property of k is that it can be used to define E as a 
Noether charge, as shown in Paper II. Quoting this result, both k and the corresponding 
energy-momentum density 

j = -T ■ k b (5.6) 

are conserved: 

V • k = (5.7a) 

V • j = 0. (5.76) 

Therefore the Gauss theorem yields conserved charges which turn out to be just the areal 
volume and energy: 

V = - / *u b ■ k (5.8a) 

E = - I *u b ■ j = [ *u ■ T b • k (5.86) 
Jt, Jt, 

where * denotes the volume form and u the unit future normal vector of an arbitrary spatial 
hypersurface E with regular centre. Paper II also showed how the Hajicek [17] energy E of 
a spherical shell of test particles of rest mass m and velocity (unit future-temporal vector) 
u may be defined analogously by 



E= I *w-T b -k (5.9) 



where the energy tensor T of the shell is 

f = m5u®u (5.10) 
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and 5 is the Dirac distribution with support on the shell. Then 

E = -mu b -k. (5.11) 

This has the same form as the standard definition of energy —mu" ■ £ for test particles in 
the stationary case [14]. 

This allows a definition of ergoregion as a region where test particles may have positive 
or negative energy E at each point, depending on their velocity u [18]. The overall sign of 
E is conventional, as in the stationary case, with E > for future pointing k in untrapped 
regions. Ergoregions are physically important because they allow energy extraction by 
the Penrose process [14] , whereby discarding negative-energy particles yields a net gain in 
energy. It follows from (5.11) that ergoregions are regions where k is spatial. Again this 
is analogous to the stationary case. In spherical symmetry, ergoregions therefore coincide 
with trapped regions. 

In summary, the Kodama vector generates a preferred flow of time in spherically sym- 
metric space-times, analogous to the Killing vector of stationary space-times and sharing 
its relevant physical properties. Moreover, it generates active gravitational energy as a 
conserved charge, suggesting the name dynamic time. This further suggests that dynamic 
gravitational entropy might be defined by analogy with the method of Wald [15]. 



VI. Dynamic surface gravity 

The above properties suggest defining dynamic surface gravity by replacing the Killing 
horizon with the trapping horizon and the Killing vector with the Kodama vector. To 
see whether this is possible, note firstly that there appears to be an ambiguity in such an 
analogy, since the Killing equation (4.1) allows the definition of surface gravity (4.2) to be 
rewritten using different linear combinations of £ • (V A £ b ) and £ • (V <8> £ b ). So one needs 
to calculate both 

k ■ (V A k b ) = (E/r 2 - Aixrw) Vr (6.1a) 
k • (V <g) k b ) = 4yrr^ b . (6.16) 

The term proportional to Vr is what is needed, since Vr = ±k b on a trapping horizon 
d±r = 0, but the term proportional to ip is generally not parallel to k. This resolves the 
ambiguity, leading to the definition of dynamic surface gravity 

E , , 

k = — — Anrw. (6.2) 

When w = 0, this has the same form as Newtonian surface gravity, or indeed Newtonian 
gravitational acceleration anywhere, with E replacing the Newtonian mass. Explicitly, the 
definition satisfies 

k ■ (V A k b ) = K Vr (6.3) 

and therefore 

k ■ (V A k b ) = ±nk b on a trapping horizon d±r = 0. (6.4) 
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This has the same form as the standard stationary definition of surface gravity (4.2). The 
new definition therefore recovers the Reissner-Nordstrom surface gravity, but generally 
differs from the definitions of both Paper I (evaluated in spherical symmetry) and Fodor 
et al. [19]. Based on the stationary case, one might conjecture that a dynamic spherically 
symmetric quantum field has a local Hawking temperature fiK/2ir, but that remains an 
open question. 

The dynamic surface gravity may equivalently be expressed as 

k = Tjdivgradr (6-5) 

where the divergence and gradient refer to the two-dimensional space normal to the spheres 
of symmetry. This may be regarded as the more fundamental, purely geometrical definition, 
with the form (6.2) following from the Einstein equations. Then outer, degenerate and 
inner trapping horizons, as defined in Section II or in Papers I or II, respectively have 
k > 0, k = and k < 0. The overall sign is conventional, but this confirms the desired 
property that surface gravity should vanish for degenerate black holes. 

The dynamic surface gravity satisfies certain inequalities involving the area and energy, 
assuming the dominant energy condition. Directly from the definition (6.2), 

Ak < 4nE. (6.6) 

This holds anywhere in the space-time, though the interpretation of k as surface gravity 
is usually intended only on a horizon. Denoting values on the trapping horizon by a 
subscripted zero and recalling Eq = |r , 



Ko < JJL (6.7) 

or equivalently k < l/2r . This has the same form as an inequality of Visser [20] for 
stationary black holes. Thus for a black hole of given area, the surface gravity has an 
upper bound. 

Moreover, Paper II established the Penrose inequality 

\P^Aq < 4tt£ (6.8) 

where E is evaluated anywhere in the untrapped region outward from the trapping horizon. 
(Recall that the outward orientation is defined for spatial and null directions). Therefore 

A K < 4nE (6.9) 

under the same conditions. If the untrapped region is asymptotically fiat, the same in- 
equality holds for the asymptotic (Bondi and ADM) energies, which are limits of E at null 
and spatial infinity respectively. This has the same form as an inequality of Heusler [21] 
for stationary black holes. Thus the total surface gravity of a black hole provides a lower 
bound for the energy measured outside the black hole. 
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VII. First law of black-hole dynamics 

The desired first law of black-hole dynamics may now be given as 

K A 1 

E ' = ^- + W V' (7.1) 

where the prime denotes the derivative along the trapping horizon, /' = z-Wf, where z is a 
vector tangent to the trapping horizon, the normalisation of z being irrelevant. This follows 
by projecting the unified first law (2.7) along z, as follows. The non-trivial part is to show 
that Aip b -z = kA' /8n. Writing z in double-null coordinates, z = z + d + + z~d-, and taking 
the horizon to be given by d+r = 0, one has = (d+r)' = z + d + d+r + z~d-d+r, which 

expands to = (d + r)' = kz~ —4tctT ++ z + , using the Einstein equations (2.10) and the 

definition of k (6.2). Therefore Ai\)° ■ z = Ag^ T ++ z + d-r = rnz~d-r = rnr' = kA' /8n. 

The term kA' /8n is exactly as expected for a first law of black-hole dynamics, with 
the derivative along the horizon replacing the state-space differential of the first law of 
black-hole statics (4.5). To complete the desired physical interpretation, the last term wV 
should be the work done along the horizon. This can be seen for the electromagnetic field, 
for which the first law becomes 

as follows from (3.3). The term S 2 V'/8tt is the electromagnetic work done along the 
horizon, in agreement with the stationary (4.5) and special-relativistic (3.5) cases. 

VIII. Quasi-local first law of relativistic thermodynamics 

Relativistic thermodynamics, as described for instance in Paper III, involves a conserved 
material current pu, where p is the density and u the unit velocity vector, 

wV = -l (8.1) 

with the physical interpretation as the flow vector of a fluid, or as giving the centre-of- 
momentum frame of a solid. The various components of the energy tensor T of the matter 
define the thermal energy density 

e = u-T^-u-p (8.2) 

the thermal flux 

Q = -±(T-w b ) (8.3) 

and the thermal stress 

t = ±T (8.4) 

where _L denotes projection by the spatial metric g + u" <E> it b , which is orthogonal to the 
material flow. Therefore the energy tensor of the thermodynamic matter is 

T = (p + e)u®u + 2u®q + T. (8.5) 
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Since specific thermal energy is e/p, this is effectively the energy tensor of Eckart [4], who 
defined specific internal energy as u ■ T b • u/p plus an undetermined constant. 

The vorticity -L(V A u") vanishes in spherical symmetry, so there exist hypersurfaces 
orthogonal to the flow vector u, labelled by a time coordinate t. Take another coordinate 
x on these hypersurfaces, orthogonal to the spheres of symmetry. Choosing vanishing shift 
vector, the line-element (2.8) becomes 

ds 2 = r 2 d0 2 + e 2C dx 2 - e 2 ^dt 2 (8.6) 

where (r, 0, () are functions of (t, x). These are comoving coordinates adapted to the fluid. 
The volume form of the hypersurfaces is given by 

where E is a region of one of the hypersurfaces and 

/ 2E\ 1 ^ 2 
a = e~ c d x r = ± I 1 + r 2 - — j (8.8) 

where the dot denotes the material (or comoving) derivative, / = u ■ V/, and the sign is 
that of d x r. The second expression gives the definition of E in comoving coordinates [8]. 
The red-shift factor a tends to one in the Newtonian limit. 

The unified first law (2.7), written in comoving coordinates, is 

d t E = Ae~ 2c: (T tx d x r - T xx d t r) (8.9a) 
d x E = Ae~ 2(f) (T tt d x r - T tx d t r). (8.96) 

For the thermodynamic matter, the relevant energy components are 

T tt = e 2 ^(p + e) (8.10a) 
T tx = -e^q (8.106) 
T xx = e 2C p (8.10c) 

where p = t£ is the radial pressure and q is the strength of the thermal flux q° = qS7x/\ Vx|. 
The heat supply as defined in Paper III is 

Q = - f*q (8.11) 
J± 

where E is a region of a hypersurface generated by flowlines of u, with volume form given 
by 

J *f = J f A ^t. (8.12) 
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Thus 

AQ = -e^Aq (8.13) 

where Af = d t f. Equivalently in spherical symmetry, Q = —Aq. Then the component 
(8.9a) of the unified first law along the flow reads 

AE = aAQ - pAV. (8.14) 

This is the desired quasi-local first law of relativistic thermodynamics. Since the quantities 
involved are local in spherical symmetry, this first law may also be written in terms of the 
material derivative as 

E = aQ-pV. (8.15) 

The first term is a relativistic modification involving a of the heat-supply term Q in the 
classical first law of thermodynamics. The last term has the same form as the hydrody- 
namic work term of the classical case, with the relevant radial pressure p and with the areal 
volume V = J s *a replacing the proper volume J E *1, again differing by the relativistic 
factor a. Recalling the role of E in black-hole dynamics, this establishes that the analogy 
between the two first laws is more than an analogy: the two first laws are different pro- 
jections of the unified first law. Both express derivatives of the active gravitational energy 
according to the Einstein equation. 

Unlike the classical first law of thermodynamics, which is an assumed principle, the 
quasi-local first law has been derived from the Einstein equations and the material model. 
Thus heat is included in general relativity as a form of energy, the quasi-local first law being 
a consequence of the relativistic definition of energy. Comparing with the non-relativistic 
first law, this indicates that the active gravitational energy E must incorporate the thermal 
energy, thereby linking thermodynamics and gravity. This may be seen explicitly as follows. 
The other component (8.9b) of the unified first law reads 

d x E=(p + s + fq/a)d x V (8.16) 

which integrates to 

E= [ *{a{p + e) + rq). (8.17) 

The integrand provides the definition of active gravitational energy density. In this form, 
E manifestly involves both thermal energy and thermal flux, though the latter disappears 
in the Newtonian limit, as follows. 

The Newtonian limit may be described as in Paper II or III by inserting factors of the 
speed c of light by the formal replacements 

(*, r) i— > (*, r) f h c _1 r p i— > c~ 2 p (£, t, E) i— > c _4 (e,r, E) q\-^c~ 5 q 

(8.18) 

then taking the limit c 1 — > 0. Then 

E = Mc 2 + H + K + U + 0(c~ 2 ) (8.19) 
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where the mass M, thermal energy or heat H, kinetic energy K and gravitational potential 
energy U are defined by 

M = 
H = 
K = 
U = 



/ *p 


(8.20a) 






/ *e 
is 


(8.206) 




(8.20c) 




(8.20d) 



All four quantities have the same form as in Newtonian theory. Thus E incorporates all 
the energies and mass present in Newtonian theory. Again this confirms the correctness of 
£asa definition of active gravitational energy. 

IX. Zeroth laws 

The zeroth law of black-hole statics states that for a stationary black hole, the surface 
gravity k is constant [13,14], a trivial property in spherical symmetry. Since a stationary 
space-time may be taken as the definition of gravitational equilibrium, this is analogous to 
the classical zeroth law of thermodynamics, which states that the temperature is constant 
in thermal equilibrium. In the Gibbs theory of thermostatics, often called equilibrium 
thermodynamics or reversible thermodynamics, this is expressed in terms of the inverse 
temperature (3 = dS/ dH. 

To see whether there is a genuine connection between these zeroth laws, restrict from 
thermodynamics to thermostatics, in which case the second law of thermodynamics 

AS > AS (9.1) 

becomes an equality. Here S is the entropy and S the entropy supply, as defined in 
Paper III. The classic relation between temperature entropy flux and thermal flux yields 

AQ = HAS (9.2) 

in spherical symmetry. Alternatively, one may regard this as defining absolute temperature 
$, as in Gibbsian thermostatics. Then the quasi- local first law (8.14) implies a quasi- local 
Gibbs equation* 

AE = a$AS - pAV. (9.3) 

If S is a function of V and E, the time derivatives may be replaced with partial derivatives 
and thermal equilibrium may be defined by constancy of 



* Authors less familiar with thermodynamics often confuse the Gibbs equation, which 
involves temperature and entropy, with the first law of thermodynamics, which involves 
neither. The two are not equivalent except in thermostatics. 
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Then 

/T 1 = a-d. (9.5) 

Therefore temperature $ is not uniform in thermal equilibrium in a gravitational field. This 
was originally predicted by Tolman [6] , who argued on different grounds that gravitational 
red-shift leads to a non-uniform temperature. In the stationary case, r vanishes and a 
reduces to (1 — 2E/r) 1 ^ 2 , which coincides remarkably with the Tolman red-shift factor. 
The above condition is more general since it allows thermal equilibrium even in a non- 
stationary space-time. Such a possibility is also predicted by relativistic kinetic theory of 
gases [22]. 

The effective temperature ad also agrees with that used for stationary black holes by 
Carter [13] and references therein. Moreover, for the case of Hawking radiation for the 
Schwarzschild black hole, the Hawking temperature 

= £ (9.6) 

is also an effective temperature, differing by the same factor from the temperature d = 
(1 — 2E/r)~ 1 / 2 f3~ 1 measured by constant-radius detectors [5]. Here units are such that 
the Boltzmann constant is unity. This confirms that active gravitational energy plays the 
role that Gibbsian thermostatics assigns to thermal energy, strengthening the link between 
gravity and thermodynamics. 

The factor k(3 = 27r/h is not determined by the above argument alone. However, it 
does show that if gravithermal equilibrium is defined by a stationary space-time with the 
above condition for thermal equilibrium, then both k and (5 are constant, which may be 
interpreted as a weak type of unified zeroth law. 



X. Conclusion 

In summary, the unified first law is just the equation (2.7) expressing the gradient of the 
active gravitational energy according to the Einstein equation, divided into energy-supply 
and work terms. Projecting this equation along the trapping horizon and along the flow 
of thermodynamic matter yield, respectively, first laws of black-hole dynamics (7.1) and 
relativistic thermodynamics (8.15), reproduced here: 

E = aQ- pV. 



In the thermodynamic case, this generalises from classical thermodynamics to general 
relativity. This shows that active gravitational energy incorporates thermal energy, as 
seen explicitly in the Newtonian limit in (8.19), thereby revealing a physical connection 
between gravity and thermodynamics. 

In the black-hole case, this shifts the paradigm from statics to dynamics, that is, from 
stationary to non-stationary black holes. The standard first law of black-hole statics, in- 
volving a state-space differential, has been replaced by a first law of black-hole dynamics, 
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involving a horizon derivative. This is ironically analogous to the paradigm shift from ther- 
mostatics to thermodynamics that has occured this century, replacing state-space differ- 
entials with material derivatives and formulating thermodynamic field equations including 
local first and second laws. The reader should be warned that this shift to genuine ther- 
modynamics has been ignored by most authors attempting to draw analogies with black 
holes, as well as by most authors of introductory thermodynamics texts. More enlightened 
references may be found in Paper III. 

Put another way, the paradigm shift is from Killing horizons and event horizons to 
trapping horizons. For instance, Hawking's second law for event horizons is replaced by 
the second law of black-hole dynamics of Paper I, which states that the area element of 
a future outer trapping horizon is non-decreasing, assuming the null energy condition. In 
spherical symmetry this reduces to 

A' > 0. 

Note that the new first and second laws of black-hole dynamics, unlike the textbook ver- 
sions, involve the same derivative. 

Gibbsian definitions of thermal equilibrium and gravithermal equilibrium have been 
given, yielding a connection between the zeroth laws, expressing constancy of surface 
gravity and an effective temperature. In contrast, the conventional second law of relativistic 
thermodynamics as described in Paper III, that the divergence of the entropy current be 
non-negative, is unrelated to the above second law of black-hole dynamics. Similarly, 
unlike the case of a quantum field on a stationary background, there will be no general 
agreement between the temperature and surface gravity, since one would expect to describe 
gravitational collapse of matter with arbitrary temperature to black holes with unrelated 
surface gravity. Moreover, it is unclear whether the gravitational entropy should agree 
with the for non-stationary black holes. 

The analogies between black-hole dynamics and thermodynamics have been known 
to be genuine in only one respect, namely Hawking's connection between surface gravity 
and effective temperature for quantum fields in the stationary case. This has remained 
so despite much effort to relate the area of a black hole to some plausible definition of 
gravitational entropy in general relativity, or indeed, in quite different theories of gravity. 
(However, see Ashtekar et al. [23] and references therein). The unified first law has now 
provided another genuine connection, neither restricted to the stationary case nor depen- 
dent on quantum theory. With hindsight, this may not seem so surprising, since both 
textbook first laws effectively express conservation of energy. However, this was just an 
analogy in the absence of any obvious link between the internal energy of classical thermo- 
dynamics and the total mass of a black-hole space-time, and indeed the latter turns out to 
be not quite right. The missing link has now been revealed as active gravitational energy. 

This provides further evidence for the physical importance of gravitational energy. In 
particular, generalising the results of this paper beyond spherical symmetry would require 
a more general definition of active gravitational energy in general relativity, on which there 
is currently no consensus. Since energy is also fundamental in quantum theory, one would 
expect gravitational energy to be important for quantum gravity. Moreover, since energy 
and time share a quantum-mechanical duality, as expressed by the Heisenberg uncertainty 
principle, this is relevant to the oft-discussed problem of time in quantum gravity. This 
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problem is neatly soluble in spherical symmetry, where the Kodama vector k encodes a 
dynamic time, related to E. This suggests that 



E = ink ■ V 

should be the energy operator of spherically symmetric quantum gravity [18]. Indeed, it 
may be argued that the correspondence principle requires this choice of time. Along with 
an area operator such as that provided by the loop theory [23], a basis for spherically 
symmetric quantum gravity may thereby be laid. 
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